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Q\ i Thin films of diblock copolymer subject to gravitational field are simulated by means of a cell 

dynamical system model. The difference in density of the two sides of the molecule and the presence 
' of the field causes the formation of lamellar patterns with orientation parallel to the confining 

Oh. walls even when they are neutral. The concentration profile of those films is analyzed in the weak 

■ segregation regime and a functional form for the profile is proposed. 

; 61.41.+e,64.60.Cn,64.75.+g 



CZ3 



I. INTRODUCTION 



Recently, there has been a great deal of interest in problems involving periodic patterns in the tens of nanome- 
ters scale, for example, light conduction by photonic crystals [Q , J osephson junctions arrays formed by granular 
. ■ superconducting materials and lithographic masks for special design chips || among others. 

In special, nanolithography has brought considerable attention to thin films of microphase separated diblock copoly- 
mers (DBCP) as they naturally self organize in periodic structures on that length scale 

The basic technique for the fabrication of those templates is the creation of a well-ordered DBCP film between two 
flat surfaces and then the transfer of the microdomains to a substrate where one of the components is removed leaving 
behind an ordered array of stripes or dots in either low or high relief. The desirable ordering, in this case, is one that 
creates a pattern on the substrate; for example, for even DBCP molecules one must have lamellae perpendicular to 
i^i - the hard boundaries. 

For that matter it is important to understand the pattern formation in confined films of DBCP since it involves 
t-H , problems not present in bulk systems. This issue has been addressed both theoretically [[|-|9) and experimentally 
Jl0p[l2| and the basic conclusions are that, when the confining walls are neutral, the equilibrium pattern corresponds 
to lamellae perpendicular to the walls and, when the substrate prefers one kind of monomer, the pattern may consist 
of lamellae parallel or perpendicular to the walls, depending on the relation between the film thickness and the bulk 
lamellar width. The later effect appears because the finiteness of the system brings about frustration and one has 
f^S ' to take into account the amount of compression or stretching of the molecules in order to accommodate a certain 
Q\ number of lamellae between the two rigid walls. 

ON ■ One important issue that has not been emphasized in the above studies is the possibility of density mismatch 
between the two parts of the molecules. At the bottom wall, as the denser part of the molecule sinks, lamellae parallel 
to the walls will form, even if the walls are neutral |L3|,|l4|. Being a bulk effect, the interaction with the gravitational 
field is capable of dramatically changing the microphase separation, even for infinite systems, as the lamellae tend 
to be aligned with the field far from the boundaries ]l3| ]. In finite systems, the lamellae perpendicular to the field 
present more diffuse interfaces and the gravitational field may completely destroy the microphase separation ]T^ |. In 
the present work we consider this problem on two-dimensional films of even DBCP molecules as we analyze the effects 
of the degree of polymerization and film thickness on frustration by means of a cell dynamical system(CDS) model. 
We simulate films both in the weak and strong segregation regimes ((WSR)and (SSR)). For the WSR we empirically 
• i-h find the one-dimensional concentration profile and study the distortion of each layer within the lamellae. 

In Sec. II, we define the model and outline the numerical scheme. Results for neutral and interacting walls in the 
presence and absence of the gravitational field are discussed in Sec. III. The effects of frustration are analyzed as they 
affect not only the size and number of lamellae, but also their internal structure. In section V, the main conclusions 
are summarized. 



II. COMPUTATIONAL MODEL 



Block copolymers are linear-chain molecules consisting of two subchains A and B grafted covalently to each other. 
Below some critical temperature T c these two blocks tend to separate, but due to the covalent bond, they can segregate 
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at best locally to form periodic structures W. Here we consider only the case of even molecules corresponding to 
lamellar equilibrium patterns. CDS models have been successfully used in several problems of phase separation 

] , so we prefer that method for the simulations. 
As usual, in this kind of description we assign a scalar variable ip(n, t) to each lattice site corresponding to the coarse- 
grained order parameter in the n-th cell at time t (time here is defined as the number of iterations). This order 
parameter represents the difference ipA — ipBi where tpA^B) is the local number density of A(B). The ingredients 
for the time evolution of ip are: local dynamics dictated by a function with two symmetric hyperbolic attractive 
fixed points, diffusive coupling with neighbors, stabilization of the homogeneous solution and conservation of ip. For 
the present problem, we also add the interaction with the gravitational field and with the confining walls. The 
conservation, when an external field is present, must be imposed by considering the Kawasaki exchange dynamics 
explicitly The detailed explanation of this model is found in for spinodal decomposition. With this, we come to 
final equation for a melt of even DBCP molecules: 

iKn, t + 1) = (1 - e)V>(n, t) + ((C (n, j; sgn[I(n, t) - t)}) [I(n, t) - t)]» , (1) 

where 

I(n, t) ee At&nh (ip(n, t)) - ip(n, t) + D [{{ip(n, t))) - tp(n, t)] + hn z + V s (n) (2) 

is essentially the chemical potential. ((*)) is the isotropic space average of *, A is a measure of the quench depth, 
and D is the diffusion coefficient. The parameter e > appears in this model to stabilize the solution ip = in 
the bulk, for e = we have a model for spinodal decomposition, in which the domains can grow without bound. 
Scaling arguments have proved that e ~ iV -2 , where N is the polymerization index f2^| , h is the gravitational field, 
which we assume is in the z direction, and n z is the z component of n. For molecules with matched densities we just 
take h = 0. A possible interaction with the walls appears via the surface term V s (n). C is the collision coefficient 
given by: C(i,j;a) = [ip c + mjj(j)][ip c — aip(i)]/ip^, where ±i/> c are the fixed points of Ataxihip — ip for A > 1 . For 
all the simulations we used the values A = 1.2 and D = 0.5, and uniformly distributed random initial conditions. 
The gravitational field, when present, is parallel to the smaller dimension. The direction normal to the field will 
be called the x direction. We consider systems with periodic boundary conditions in the x direction and hard walls 
in the z direction, separated by a distance L z . At the hard walls we impose no flux boundary conditions in the 
form: [1(2 + 1) - I (z)] boundaries = 0. 



dynamics due to their computational efficiency and versatility [TL5 



III. RESULTS 



In order to understand the effect of confinement on the lamellae width, we must first determine its bulk value. For 
that matter, we ran simulations on 512 x 512 lattices with periodic boundary conditions for different values of e, and 
h = 0. The resulting isotropically striped pattern was then Fourier transformed, and the bulk lamellar width Wb was 
measured in a standard way. Defining one lamella as ABBA we have: 

where 

_ JS(k,ac)kdk 

{)eq ~ fS(k,<x>) ■ [) 

S(k,oo) is the circular average of the structure function S(k,t) — \ip(k, t)| 2 , calculated at large times, that is, when 
the value of (k) approaches a constant value. 

Just to make sure that the bulk lamellar width as measured above was not affected by the interface bending of 
the disordered pattern, we also measured Wb in 32 x 128 systems with hard neutral walls and zero field, or matched 
densities. As expected the equilibrium configuration correspond to lamellae normal to the hard walls ||[^| as in Fig. 
|l|(a) . Wb was then measured using the one-dimensional structure function for each line and finally averaging along the 
z direction. The values of Wb found in both determinations agree, so we conclude that the excessive interface curving 
of the disordered pattern does not affect the lamellar width. Since the disordered patterns are easier to obtain, we 
will consider the lamellar width obtained from them as our bulk equilibrium value Wb- 
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The results below correspond to simulations with V s = (neutral walls), h ^ (mismatched densities) and h = 
(matched densities), and V s ^ (interacting walls) and h = 0. As will be seen, different patterns regarding the 
lamellae orientation appear: lamellae normal or parallel to the hard walls and a mixture of both. 
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( Q )h=0;e = 0.004; V s = 



^^^^ 

(b)h = 0.01;e = 0.0014; V s = 



(c)h=0.01;e = 0.01; V s = 



(d) h=0; e = 0.01 ;V, = 0.01 

FIG. 1. Equilibrium patterns for confined films with L z = 21 and L x = 256 (only the first 150 columns are shown), (a) 
Neutral walls, and matched densities. The lamellae are normal to the walls with the bulk periodicity, (b) Neutral walls and 
density mismatch, e in the SSR. 1.5 lamellae are accommodated parallel to the walls, (c) Neutral walls and density mismatch, 
e in the WSR. 2.5 lamellae are formed parallel to the walls. The lamellar width is 8.401, smaller then the bulk value of 9.422. 
(d)Interacting walls and matched densities, e in the WSR. The lamellae are also parallel to the walls but are more segregated 
than in (c). 



A. Neutral walls 



We focus now on hims with mismatched densities (h ^ 0), confined by neutral walls (V s =0). In this situation we 
observe patterns of lamellae parallel to the hard walls, or a mixture of wetting layers on the hard walls and lamellae 
normal to the walls in the center part of the film. First we analyze the case of lamellae parallel to the walls only. 
Due to the density difference of chains A and B, the denser part (say A) will be at the bottom, and the less dense 
part (say B), at the top. For a blend of two homopolymers A and B, the film would have the lower half filled with 
A and the upper one with B. The covalent bond between A and B parts hinders this complete separation and forces 
the alternation of A-rich and B-rich microregions that will then have thicker interfaces due to the interpenetration of 
domains ]l4|] . In the extreme case, the existence of a density mismatch may completely destroy the segregation of A 
and B. The number of alternating lamellae will depend on both chain size, e, (Fig. [lj) and the separation between 
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walls, L z , (see Figure g). Also, the equilibrium patterns will always have m + 1/2 lamellae, where m=0,l,2. . . . 

P I 



(a) L z = 29 

Ji AIU>~ 

(b) L z =26 



(c) U =20 




(d)L 2 =17 



(e) L z =14 




(f) L z =9 



(d) L z =5 

FIG. 2. Equilibrium patterns for films with h = 0.01, e = O.Ol(WSR) and variable width. L x = 256 but only the first 150 
columns are shown. As the width L z is decreased, the him goes discontinuously from m = 3 to m = lamellar patterns, (b), 
(d) and (f) show transition patterns with mixed parallel and perpendicular lamellae. 

As we vary L z for a fixed e we clearly see the effects of frustration. Figure || shows the transitions from m = 0, 
to m = 1, m = 2 and m = 3 patterns as L z is changed from 5 to 29 for e = 0.01 and /i = 0.01. The transition 
patterns are frustrated and present lamellae normal to the walls in the central region. Since full lamellar patterns are 
essentially one dimensional, we define the average concentration profile, {ip{n z )) Xl as the average over the x direction 
of the vertical variation of ip. Figure || shows the behavior of (i[>(n z )) x , for three different situations, for now we are 
interested in cases (a) and (b) only. Figure ||(a) corresponds to the profile for e = 0.01, h = 0.01, V s — and L z = 21. 
If we try to fit a sine function to that profile, we see that the fitting will miss only the wetting layers. From this we 
conclude that the system is in the WSR so that the inner layers can be described by just one Fourier component p3| . 
The wetting layers have an enhanced concentration due to gravitational field and the presence of the wall: the bottom 
and top AB layers are considerably stretched by the effect of buoyancy and do not experiment the penetration of 
other layers, resulting in an excess of A at the bottom and of B at the top. 
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(a)h = 0.01;e=0.01; o= 



5 10 15 20 

(b)h = 0.01; e= 0.0014; a=0 



5 10 15 20 

(c) h = 0; £=0.0014; o= 0.01 



FIG. 3. Average concentration profiles for a film with L z = 21 and L x = 256. The continuous line correspond to a fitting 
using Eq. (^) (a)Neutral walls, mismatched densities. The profile is well fitted by a sine function plus exponential enhancement 
at the hard walls, (b) Neutral walls, mismatched densities and molecules larger than in (a). From the fitting it is clear that 
more than one Fourier component must be considered, indicating that the film is already in the SSR; (c) Interacting walls 
(V s (l) = — o = —Vs(L z )) and matched densities. Although similar to the profile (a), we notice that Eq. (0) is not adequate to 
describe this concentration profile. 

A correction for this effect led us to the tentative function 

ip{x) = (-l) m+1 ? ? sin g iE + 2Ce- /3i sinh2/3:z; for -L/2 <x< L/2, (5) 

which fits very well the profile in Figure |](a). For films in the WSR we found that the fitted value for q is indistin- 
guishable from 27r(m + 1/2)/ L z , so we define the average lamellar width W, directly from the fitting, as 2n/q. As 
will be seen below, the interaction with the gravitational field causes a distortion within the lamellae regarding the 
width of the A and £?-rich layers, but, if considered as a unit, all the lamellae have a width very close to the average 
value. The transitions between consecutive values of m as we vary L z is shown in Figure ||. From this figure we see 
that discontinuous transitions occur from a pattern in which the lamellae are stretched, compared to its bulk state, 
to a compressed state with one more lamella, as L z is increased. The regions between steps of fixed m correspond to 
transition patterns in which lamellae normal to the walls form in the center portion of the film. 
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FIG. 4. Total number of lamellae, n, as a function of the film thickness L z for a film with mismatched densities, confined by 
neutral walls. The solid points on the steps represent lamellar patterns with n lamellae parallel to the hard walls and between 
consecutive steps the film is in a mixed configuration with horizontal and vertical lamellae. The vertical dashed lines correspond 
to the film thickness adequate to accommodate the corresponding number of lamellae but with the bulk width Wb- We see 
that as the film width increases, discontinuous transitions between n stretched lamellae and n + 1 compressed lamellae occur. 
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We can, alternatively, fix L z and vary e, which corresponds to fixing the film width and varying the bulk lamellae 
width. For 0.006 < e < 0.018 we obtain patterns with m = 2, in the range 0.004 < e < 0.006 again we observe a 
transition pattern, and decreasing e further we find that a m = 1 pattern appears. Figure || shows patterns with 
m=l, to=2 and in the transition region. The analysis of the transitions in this case is more complicated since for 
e < 0.004 the system is no longer in the WSR as can be seen from the fitting of Eq. (^) in Fig. ||(b). It is clear that 
other Fourier components need to be included in this case. 



= 



0.01 



£ = 0.005 



e = 0.0014 

FIG. 5. Equilibrium patterns for confined films with L z = 21 and L x = 256 (only the first 150 columns are shown) three 
different values of e. As in the case of variable film width, the number of lamellae varies discontinuously and transition patterns 
with mixed orientation lamellae appear. For e = 0.01 2.5 weakly segregated lamellae are formed, the concentration profile can 
be well fitted by Eq. |^. e = 0.005 produces a transition pattern and e = 0.0014, a more segregated pattern with 1.5. 

The accommodation of the lamellae distorts their widths non uniformly as may be easily checked from the plot of 
the width of each individual A-rich and B-rich layer. Figure || shows the behavior, as a function of e, of the widths 
of the first A-rich layer that wets the bottom wall the first _B-rich layer connected to the bottom wetting layer 

(1V2) and one quart of the central lamella (w c — W/A). For the sake of comparison, the variation of the bulk width, 
Wb = Wb/i, of A-rich (or B-rich) layers is also plotted. Although the variation is small compared to the bulk behavior, 
we see that W2 < w c < w\ consistently This happens because there is a reduction in the number of A and B contacts 
in the first layer (for the lack of neighboring molecules from below) and an increase in the next one because the 
gravitational field shifts the A parts downwards. As we separate the regions where the lamellae are compressed and 
stretched as compared to the bulk, we notice two different situations. In the compressed region, wi increases as N 
increases (e decreases), due to a compression of the internal layers caused by the greater stretching of the surface 
layers which produces an increase in the internal pressure (for L z fixed). We expected the inverse effect to happen 
when N was reduced in the stretched region: the internal layers would shrink producing a tension that would stretch 
the surface (larger effect) and second layers (smaller effect). But, in fact, we observe a drastic reduction of the second 
layer acting as a tension center for surface and central layers (Fig. 0). 
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FIG. 6. Variation of ^4-rich and B-rich layers as a function of e for a films with 2.5 lamellae. Wi, W2 are the widths of the 
lowest A-rich and B-rich layers, w c is 1/4 of the central lamella and Wb is 1/4 of the bulk lamellar width. The vertical line 
indicates the separation between regions where the film is in compressed and stretched states. 

It is clear that the above effects are meaningful only for thin films. The transition from this to the bulk behavior 
may be observed as we analyze Wi, u>2 and w c as a function of the film thickness L z . If the bulk behavior prevails, 
W2 ~ wi ~ W2 ~ wi, w L z /m. As we increase L z and observe films with increasing number of lamellae, we find 
wi — > W2 — * w c — > Wb- On the other hand, the slope a of each group of w values is proportional to to -0 8 instead of 
?7i _1 , which reflects the different behavior of each layer of DBCP under stretching or compressions. 
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FIG. 7. Variation of A-rich and B-rich layers as a function of L z for a films with e = 0.01. We use here the same notation 
of Fig. ^| Each group of points corresponds to lamellar patterns with (m + 1/2) lamellae. As L z , and correspondingly m, 
increases, the film behaves more like a bulk sample in the sense that the distortion of lamellae is less significant. The slope a 
of each group is proportional to m -0 ' 8 , for larger values of L z a crossover to the bulk behavior a oc m _1 is expected 



B. Interacting walls 

The effect of surface fields in the formation of lamellar patterns has been extensively studied — . Our goal here 
is to compare the effect of the surface and bulk fields, so we consider only a film of DBCP molecules with matched 
densities confined by interacting walls in such a way that the bottom wall attracts the denser component and the top 
wall prefers the less dense component, As will seen below, in many ways this choice of walls produces a pattern is 
similar to the one obtained with neutral walls and a density mismatch, but the two situations are, in fact, different. 

The above interaction with walls may be simulated by choosing the surface interaction as: V s — a for z = 1 and 
V s = —a for z = L z . The equilibrium pattern obtained for e = 0.01, L z — 21, h = and a = 0.01 is very similar to the 
one with the same values of e and L z , a = and h = 0.01, as both present 2.5 lamellae parallel to the walls (see Fig. 
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|l| (c) and (d), respectively). The first distinction appears in the segregation of domains: it is clear that the pattern in 
Fig. H(c) is less segregated due to effect of interpenetration of domains driven by the gravitational field. As we try to 
fit the concentration profile with Eq. (|J) we notice that the patterns with density mismatch and surface interaction 
are also a little different in the surface region, so we conclude that Eq. (||) is a good fit for the concentration profile 
of lamellar patterns with density mismatch and neutral walls only. A substantial difference appears for larger values 
of h and a. Figure || shows patterns with the same value of e = 0.01 and L z — 21 but one with surface field only and 
the other with density mismatch only. In this case, the lamellar structure still exists for tr = 0.04 but here, for h = a 
the lamellar structure is completely destroyed in the center of the film. 
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(a)e = 0.01; h = 0; V s = 0.04 



\ • 4 A A 4 4 




(b) e = 0.01; h = 0.04; V s = 



FIG. 8. A comparison between bulk and surface interactions, (a) Interacting walls and matched densities. Although the 
surface interaction is stronger than the one considered so far, the only noticeable difference is the segregation of the wetting 
lamellae, (b) Neutral walls and mismatched densities. Increasing the value of the bulk field h by the same amount as the 
surface field the observed pattern changes dramatically: instead of a lamellae pattern we observe a frustrated mixed orientation 
pattern. 



IV. CONCLUSIONS 



In this paper, we study the effects of surface and bulk (gravitational) fields coupled with hard wall restrictions on 
the lamellar pattern formation of diblock copolymer systems. We find that the two are the predominant factors to 
determine the final equilibrium pattern: lamellae tend to form normal to the field and their number is determined 
by the ability of the system to resolve the frustration caused by the confinement. Unresolved patterns present a 
mixture of wetting lamellae normal to the field and lamellae parallel to the field in the central part of the film. The 
gravitational field also distorts the periodicity of the lamellar pattern. The bottom A layer is larger than it would be 
if placed in the central part of the film. On the other hand, the next B layer is narrower, in such a way that the first 
lamella, defined as the sequence ABBA, has a width very close to the central lamellae. 

We obtained a good fit for the average concentration profile in the WSR by using a trial function which consists 
of the superposition of a sinusoidal function, characteristic of the WSR, and exponential functions for the enhanced 
concentration of the wetting layers. 
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